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Abstract 

We describe a method to accelerate the numerical computation of the coef- 
ficients of the polynomials Pk{x) that appear in the conjectured asymptotics of 
the 2A;-th moment of the Riemann zeta function. We carried out our method to 
compute the moment polynomials for k < 13, and used these to experimentally 
test conjectures for the moments up to height 10^. 



1 Introduction 

For positive integer k, and any e > 0, Conrey, Farmer, Keating, Rubinstein, and 
Snaith conjectured |CFKRSj that 

|C(l/2 + tt)\''dt = r Pk (log(t/(27r))) dt + (ri/2+^) , (1) 



'0 JO 

with the constant in the O term depending on k and e. 

In the above equation Pk is the polynomial of degree k'^ given implicitly by the 
2/c-fold residue 
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G{zi, . . .,Z2k)^^{zi, ■ ■ .,Z2k) 
2k 

i=l 

X 61^^=1^'"^'+* dzi... dz2k, (2) 
with the path of integration over small circles about Zi — 0, where 
A{zi, ...,z^)= J] {^Zj - Zi) = Izf' 1^^^ 

l<i<j<m 

denotes the Vandermonde determinant, 

k k 

G{zi, Z2k) = Ak{zi, Z2k) 1111'''^-^ + ^*" ^•?'+'^)' 

i=l j=l 

and Ak is the Euler product 

Ak{Zi, . . . , Z2k) 

Here e(^) = exp(27r«^). 

We denote the coefficients of Pk{x) by Cr{k): 

Pfc(a;)=:^c,(A;)a;^'-'-. (5) 

r=0 

In order to arrive at this conjecture, CFKRS considered a more general moment 
problem with 'shifts'. Because this general setting was central to our computation, 
we describe their conjecture with shifts below. Write the functional equation of zeta 
as 

C(s)=x(s)C(l-5), 

where 

X{s) :=7r-^/2r((l-.)/2)/r(./2). 
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It is a little more convenient to work with the Hardy Z-function, whose functional 
equation and approximate functional equation are expressed more symmetrically 
than that of the zeta function. It is defined as 

Z{s) = xisr'/'cis), 

and satisfies 

i) Z{s) = Z{1 - s) (because x(s)x(l - s) = 1), 
n) Z{l/2 + it) e R for t e M, 
iii) \Z{l/2 + it)\^\C{l/2 + it)\. 
CFKRS took as their starting point the shifted moments: 

M{ai, a2k) ■■= [ Z{l/2 + it + ai)--- Z{l/2 + it + a2k) dt, 
Jo 

where aj G C are distinct and satisfy — 1/4 < Kctj. When a — the integrand is 

|C(l/2 + ^t)p^ 

Substituting the approximate functional equation into each factor of the above 
integrand 

s = a + it, < a < 1, CFKRS apphed the following heuristic steps: 

a) Ignore the 0{t~'^^'^) and expand the product to get 2^^ terms, each a product 
2k sums. 

b) Of the 2^^^ terms, only the terms with the same number of s's and 1 — s's 
contribute to the asymptotics. Reasoning: x(s) is highly oscillatory, so cancel- 
lation occurs unless each s gets paired with a 1 — s. 

c) For any such term, only the diagonal ('mim2 . . . nik — nin2 ■ ■ ■ rik) contributes 
when the sums are multiplied out. 

d) Extend the truncated diagonal sums to infinity, replacing the sums that diverge 
with their analytic continuation (the assumption we stated earlier, —1/4 < 
^aj, is used when obtaining the analytic continuation of the diagonal sums). 
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The steps in this heuristic recipe are not justifiable. In fact, the terms that are 
dropped cannot be neglected individually, and it appears that some sort of can- 
cellation takes place amongst these terms so that, in the end, the above steps do 
apparently result in a correct conjecture. 
Let 

H{zi, . . .,Z2k;x) := exp (^^Zj - Zj+k^ G{zi, Z2k)- (6) 

The first conjectured asymptotic formula of CFKRS, which we refer to as the com- 
binatorial sum, for shifted moments reads: 

M(ai,...,a2fc)~ /" Pfc(a,log^) rft, (7) 

where 

Pfc(a, x) = ^ H{a^(^i), a^{2k)] x), (8) 

and S is the set of i^^^ permutations such that cr(l) < ■ ■ ■ < (j{k) and (j{k -|- 1) < 
■ ■ ■ < a{2k). The terms in this set correspond to the number of ways to select, from 
approximate functional equation for the same number k of s's and (1 — s)'s. 

CFKRS, also expressed their sum of {^^^ terms as a 2A;-fold residue. We reproduce 
their second formula for the shifted moments: 

(-1)^ 1 / 
''^(^-^^^ k}- {2my^ f 



While the (^^^) terms of the combinatorial sum ([8]) for Pk{a,x) each have poles 
of order /c^ at a = 0, the above is analytic in a neighbourhood of a = which shows 
that these poles must cancel. Working with shifts allowed CFKRS to get beyond 
these poles. 

To get formula ([2]) for Pk{X), set a = in ([2]). Even though the formula conjec- 
tured by CFKRS is complicated, it does seem to correctly predict the moments. Some 
evidence, both theoretical and numerical, in its favour was presented in |CFKRS] 
and |CFKRS2] . For instance, the formulas predicted by CFKRS match with known 
theorems, including lower terms, for k = 1,2. CFKRS also computed the moment 
polynomials for k <7, and tested the moment conjecture for T roughly of size 10^. 



G{Zi, . . . , Z2k)^{Zl, . . . , Z2k)' 
2k 2k 

i=i j=i 

xe2 -^'^i^' '+''dzi . . .dz2k- 



(9) 
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In this paper we describe a 'cubic accelerant' variant of the second method pre- 
sented in |CFKRS2] for computing coefficients of the moment polynomials. This 
allowed us to extend tables of the Cr{k) for k < 13, and, < r < A;^. We also 
computed many of the coefficients to greater precision. Finally, we used our tables 
of coefficients to test the moment conjectures up to A; < 13 and T up to 10*^. 

Going up to A; = 13 is substantial, because computing the coefficients to D digits 
accuracy involves evaluating k"^ sums (one for each r), each sum involving 2k choose 
k terms (10,400,600 for k = 13), with working precision of roughly D x k"^ digits 
accuracy. For example, about 2000 digits are required for A; = 13 with desired 
precision of 12 digits. The process is made even more challenging, by the fact that 
each term involves a complicated infinite multivariate product over primes. The 
RAM requirements are also large because of the number of terms, (^^) , that we are 
computing/updating one prime at a time to high precision. 

2 Computing the coefficients of Pk{x) 

In |CFKRS2] , two methods were described for numerically computing the coefficients 
Cr{k) of the moment polynomials Pk{x). 

Their first method involved expanding, in a multivariate Taylor series, the inte- 
grand of ([2]), and working out a technique for expressing the resulting residue (for a 
general k), giving formulas for Cr{k). For example. 



Co (A;) 




ci{k) 

C2{k) 



Co(A;)2A;2(7A; + 5fc(l;)) 
Co(A;)A;2(A;-l)(A; + l) 

X (2(5fe(l; ) + ^kf - 7^ - 271 + 3^(1, 1; ) - 5,(1; 1 



)) 



(10) 



where 



ak = ll{l-p-'f 2Fi{k,k;l;l 



Iv) 



(11) 



V 



the 7j's 



defined by 



s C(l + s) = 1 + 70S - 71^2 + f - f + ■ ■ ■ ' 



(12) 
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and 



>^ fclog(p) \og{p)k 2Fi{k + 1, fc + 1; 2; 1/p) 
^ p-1 p2F,{k,k;l-l/p) 



p 



(13) 



, . ^ >p ( \og{pfk\F,{k + 1, + 1; 2; log(p)^(^+^)2Fi(fc + 2, + 2; 3; ' 

' p\F,{k,k-l-l/pY p^2F,{k,k-l-l/p) 

p\og{pf (\og{pfk'^2Fi{k + 1, + 1; 2; l/pf 



' ^ ^ {p~iy \ p\F,{k, k- 1; l/pf 

log{pf2Fi{k + l,k + l;l;l/py 



P2F,{k,k;l;l/p) 
^ ^ kplogjpf f \og{pfe2F,{k + l,k + l;2;l/pf 
^ (p-1)^ V p^2F^{k,k■l■l/p)^ 

log(p)^ Ct') + 2, fc + 2; 3; 1/p) log(p)2A: 2Fiik + 2, A: + 1; 2; 1/p) 

p2 2Fi(A;,A;;l;l/p) p2Fi(A;, A;; 1; 1/p) 

with 2F1 Gauss' hypergeometric function. 

These formulas quickly get much more complicated. In practice CFKRS were able 
to use this method for r < 9 and compute numerical approximations for all the coef- 
ficients of, for example, Ps^x). One advantage of these formulas, expressed as sums 
over primes, is that one can apply Mobius inversion to accelerate the convergence of 
these sums, and obtain high precision values of the coefficients. 

It soon became apparent |CFKRS] from numerical values of Cr{k) that the lead- 
ings coefficients of Pk{x), i.e. associated to the larger powers of x, are very small 
in comparison to the lower terms. Thus, in order to meaningfully test the moment 
conjecture for zeta, which involves the moment polynomial evaluated at the slowly 
growing function log(t/27r) (and this hardly changes over the range of t in which we 
can gather significant data for ({1/2 + it)), one needs many coefficients of the mo- 
ment polynomials. See also |HR] which discusses the uniform asymptotics of these 
coefficients. 

Consequently, a second practical method, relying on the combinatorial sum ([8]), 
was developed for computing numerical approximations for all k"^ coefficients of the 
moment polynomial Pk{x). 

We detail our computational approach, implementation, and numerical results in 
the next two sections. 
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3 Our numerical evaluation of Cr{k) 

The polynomial Pk{x) given by ([2]) is the special case «! = ... = = of the 
function Pk{a,x) in ([9]). CFRKS's second method for computing the coefficients of 
Pk{x) relies on their equation ([8]) for Pk{a,x). However, the terms in ([8]) have poles 
if the ctj's are not distinct, coming from the product of k"^ zetas, 

k k 

l[l[C{l + z,-z^^,), (14) 
1=1 j=i 

that appear in the function G. So we cannot simply substitute aj = 0. 

Instead we take the limit as aj while making sure that all the a^-'s are distinct. 
Because of the poles, each individual term in ([8]) becomes very large when a is small, 
and high precision is needed to see one's way through the resulting cancellation of 
the poles as we sum across the (^^) terms of the combinatorial sum. 

More precisely, consider 

(A: \ k k 

2^^j~ ^j+k I Ak{zi, Z2k) n n ^'^■^ + ~ ^^■+'^)' 
1 / i=l j=l 

(15) 

and let 

ej = je, (16) 

where e G C is a small number. In practice e was of the form 10^^ for some positive 
integer D. 

Using (IHl) we obtain 

Pk{x) = limy^i^(e^(i), . . . ,e^(2fc);a;). (17) 

As in |CFKRS2] . we expand exp in its Taylor series, and pull out the coefficient 
of x'^ to get 

^^rik) = lim^g,(e.(i), . . . ,e^(2fc)), (18) 



where 



Hr{zi, . . . , Z2k) = \ - Zj+k I Afc(zi, . . . , Z2k) JI 11 ^^^^ + ~ ^^^^'^^ ^^^^ 

1 / i=l j=l 
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Notice that, as a function of e, Hr{€cr{i), ■ ■ ■ ,€cr{2k)) has a pole at e = of order r, 
because the first factor above cancels k"^ — r of the A;^ poles of the double product 
of zetas. These poles must cancel when summed over permutations a, otherwise we 
would not obtain the Ihs, Cr{k), as e — )■ 0. Therefore, because the sum in f|T8|) is 
analytic about e = 0, we may write 




^rik) = _ 1^ 1^2^ g.(Mi), ■ ■ ■ , 6.(2.)) 1 +0(|e|) 1 , (20) 

with the implied constant in the remainder term depending on k and r. 

One complication in evaluating the above for a given k and e is that Ak{zi, . . . , Z2k) 
is expressed as an infinite product over primes as described by @. 

While CFKRS used a 'quadratic accelerant' for evaluating the multivariate Euler 
product, we implemented a cubic accelerant. This has the advantage of allowing us 
to truncate the Euler product sooner. 

To evaluate Ak{zi, . . . , Z2k), we break up the product over primes into p < P and 
p > P, where P is a large number. For the first portion p < P, we use the following 
identity, derived in Section 2.6 of [CFKRS], 

k 

A>(.........>)^nEn "^'i ■ pi) 

to numerically compute the local factor of ^^.(zi, . . . , Z2k) for specific values of p and 

Zl, ■ ■ ■ , Z2k- 

Some care must be taken to account for the fact that individual terms in this 
identity also have poles. While these poles cancel out when summed over j, see the 
paragraph following equation (2.6.16) in |CFKRSj . they cause some additional loss of 
precision in our application. We are evaluating Ak{zi, . . . , Z2k) at distinct, but small 
values of zj. Therefore, when evaluating the sum over j, additional cancellation and 
hence loss of precision occurs affecting the leading {k — 1)D digits, where e ~ 10~^, 
from the poles of order k — 1 of the individual terms summed in fl2ip . 

For the contribution of the second portion p > P, we approximate each local 
factor appearing in by a product of zeta functions that captures the terms up to 
degree three in 1/p of its multivariate Dirichlet series. A cubic approximation can be 
obtained by first substituting uj = and Wj = p-^/^+^fc+j into the local factor 

of (H, 

n (1 - u,w,) I n(l - u,e{e))-\l - w,e{-e))-' dd, (22) 

i,j=l "^0 j=l 
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and then working out the terms, in the multivariate Maclaurin series, up to degree 
six, in Uj and wj. 

Notice that the integral over 6 pulls out just the terms with the same number of 
m's and w's. This results in monomials only of even degree appearing. The integral 
of any other term, which does not have the same number of m's and w's, is zero 
because it contains a non-zero integer power of e{9). 

Also, observe that the local factor of ( l22l) is symmetric in the m's and, separately 
in the w's, meaning if the Ui's are permuted the expression remains invariant and 
similarly for the Wj's. Also, it is symmetric with u and w, i.e. if all the m's and w's 
are swapped the expression remains the same. 

Therefore, to get terms up to degree six, we can determine the coefficients of 
representative terms involving ui,U2,U3 and wi,W2,W3, and then symmetrize the 
resulting expressions over all the m's and w's. More precisely, to get all terms of 
degrees 2, 4, and 6, it is sufficient to consider only the monomials: uiWi, U1U2W1W2, 
uiU2wl, ulwl, U1U2U3W1W2W3, U1U2U3WIW2, uiU2U3wf, u\u2w'lw2, u{u2w\, and u\w\ 
instead of every possible monomial, and then exploit symmetry. 

Finally, we can simplify further. The integral over 6 simply plays the role of 
pulling out terms with the same number of m's and w's. So, instead of fl^2]) . we can 
work more directly with the function 

k k 

n (1 - n(i - u^)-\i - (23) 

i,j=l j=l 

The multivariate Maclaurin series of the above coincides with that of fl22|) for those 
terms that have the same number of m's and w^s. Furthermore, because we are 
focusing just on terms involving ui, U2, M3 and Wi, W2, W3, we can set Uj = Wj = for 
all j > 4. Finally, to get terms up to degree six with the same number of m's and 
w's, we can expand each factor in the denominator as a geometric series of degree 3. 
We therefore consider: 

3 3 

JJ (1 - UiWj) + Uj + u] + + Wj + w] + (24) 

ij=l j=l 

We expand out the above and tabulate, in Table [H the coefficients for representative 
monomials with the same number of u's and ly's, up to degree 6, in the multivariate 
Maclaurin series of the above function, and, hence, equivalently, in f l22|) . 
Therefore, symmetrizing, we have that f l22|) equals 

1 — Ui^Ui^w j{w j.^ + 4 Ui{ai.{ai.^w j{u) j^w j.^ 

I<i-|^<i2!^fe I<i-[^<i2<i3<fc 
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monomial 


ropfRripTit 

vj V_jXXX vjX V^XX U 




n 

w 


UIU2WIW2 


1 

-i 


2 

UiU2W^ 


U 




n 

u 


U1U2U3W1W2W3 


4 


U1U2U3WIW2 


1 


U1U2U3WI 





U\U2WIW2 





U\U2W\ 





3 3 

u'lw'l 






Table 1: The second column lists the coefficients that appear with representative 
monomials, up to degree 6, in the multivariate Maclaurin expansion of (122!) . 



Yl Ui^Ui^Ui^w^Wj^ + ului^Wj^Wj^Wj^ + ... (25) 



1 <i]^ <i2 <^3 l^i-l^i2—^ 



Undoing the substitution for Wj's and w/s, gives the following expansion for the 
local factors in (HI): 



1 - E 



l<Jl<j2<fe 



+ 4 5^ 



-3 Zi-^ Zi^ Zi,^+z^.^j^+z^.^j^+z^.^j^ 



I<i]^<i2<i3<fc 
l<il<j2<^3^* 



+ Y P 

l<i-^<Z2<i3<fc 



3 '2,Zi-^ Zi2 -\-Z}^j^jj^ -\-Zf^-\-j2 -\-Z]^^j^ 



1 <i]^ 7^*2 — ^ 
1<J1<J2<^3<^ 



+ . . . , (26) 
which we then approximate by the following product: 

( 1 — ""^'2 +^fe+Jl +2fc+J2 ) 



l<z-[^<i2<fc 
1<J1<J2<'= 

X 



l<jl<J2<J3^^ 
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X 



(1 — p~3-2;ij-2i2-2i3+22fc+j^+2fc+j2^-l 



l<i-|^<Z2<''3^fc 



l<i ]^ 7^*2 — ^ 



(27) 



The last step can be seen by expanding each factor in a geometric series and compar- 
ing the terms, up to those containing a with those in f l26|) . We also remark that, 
had we wanted a quartic approximation, then slightly more care would be needed as 
the first product above would, on expanding in geometric series, interact with the 
quartic terms. 

Note that, when all Zj = 0, the above gives an approximation for the local factor 
with a remainder term of size Ok{p~^)- Hence, the accumulated remainder for p > P 
in this approximation, when Zj =0, is 

+ Oip-')) « J2p^' « ^ogiP)-'p-', (28) 
p>p p>p 

with the implied constant depending on k. 

The product in ( j2711 allows us to approximate the tail, i.e. for p > P, of 
Ak{zi, . . . , Zk) in terms of zeta: 



p>PiJ=l j=l 



n l<n<i2<k C(2 + ^ii + — Zk+j^ — Zk+j2) ^ 
l<jl<j2<fc 

Up<P n i<n<'2</^ (1 - p-^-^n--^2+-k+n+->^+n ) 

l<.7i<i9<fc 



<ii<i2<'' 

4 



Y\ i<ii<i2<'3<'= C(3 + ^ii + ~^ ^i'i ^k+ji ^A;+j2 ^k+j: 
^ I<jl<j2<j3<'= 

n <pni<n<-2<«3<* (1 -p-3^^n-^«2-^«3+^fe+n+^fc+.2+^'=+.3)-4 

'- ^<jl<32<ri<k 

I<!l<i2<i3<fc C(3 + Zi^ + Zjj + 2:43 — 2Zk+ji ~ ^k+ji) 

n <pni<n<-2<«3<fc(l -J5-^~''i-"'2-^«3+2-fe+n+^*+.2)-l ^ ' 

n i<ii7^i2<* C(3 + 22j^ + 2;j2 — Zk-\-ji — Zk+j2 — Zk+jg) 

^ i<ii<i2<i3<fc 

n„<pn l<n7^«2<fe (1 -p'3-22,^-2,2+2fc+,^+2fc+J2+2fc+J3^_l• 
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4 Implementation and tables of coefficients c^(/c) 

Our code was implemented in C++ using the GNU MPFR library jFHLPZj . along 
with Jon Wilkening's C++ wrapper for MPFR. MPFR is based on GMP, the GNU 
multiprecision library. We used gcc, the GNU C compiler, to compile our code, with 
the '-fopenmp' option in order to enable the use of OpenMP directives in our code. 
This allowed us to carry out some of the key steps in parallel for a given k, using 
several cores of our machine. 

For each k, we selected a precision, specified by the number of digits desired, 
'Digits', for the final output, and let e = 10"^^^'*^. For example, we used Digits= 25, 
i.e. e = 10~^^, for k = 4. We then put ej = je, I < j < k, and set about computing 
the sum (120]) . using our cubic multivariate approximation for the tail of the infinite 
product, i.e. equation fl2^ . 

Observe, in (fT9l) . that the dependence on r manifests only at the factor: 



Therefore, we were able to store and recycle all the other quantities across < r < k^. 

We record one important hack that we used several times in our program. While 
the double product of zetas in (IT^ involves k'^ factors, many of these are repeated 
since there are just 4A; — 2 possible values of — = (a — b)e, where a, b are distinct 
integers in [1, 2k]. The same holds, for each p, in the double products in fl2T]) . 

Likewise, while the products in fl29|) involve up to 0{k^) factors, these appear 
with multiplicity, and there are just 0{k) distinct factors. This is true both for the 
product of zetas in the numerator, and also, for each p, the factors that appear in 
the denominator. 

We exploited these multiplicities by computing and storing a table of the distinct 
values of zeta that appear and, for each p, of the distinct factors that occur. Further- 
more, we took advantage of the fact that the powers of p that occur, other than 1/p, 
are of the form p™^, where m G Z, and thus computed and stored them by repeated 
multiplication of p*^ and of l/p*^. 

To account for the high amount of cancellation that occurs as a consequence of 
the poles of the individual terms in fl20l) . we let our working precision be equal to 



and carried out our computations using these many digits. The k'^ was to account 
for the largest order poles occurring in H^, when r = /c^, of order k"^. While we could 




(30) 



WorkingDigits = (k"^ + k — 1 + 6)Digits, 



(31) 
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have gotten away with less precision for smaller r, we recycled most of the computed 
quantities across all r. The k — 1 acounts for cancellation amongst the poles of the 
terms in the sum over j in (|2T1) . Finally we needed to have some working precision 
left over, after all the cancellation, to capture Cr{k) to Digits precision. The +6Digits 
was chosen to give us some leeway. For example, we had WorkingDigits= 625 for 
k = A and Digits= 25, and WorkingDigits= 2244 for A; = 13 and Digits= 12. 

Note that using a specific desired precision does not necessarily result in that 
precision being achieved as one also needs to take into account the implied constant 
in the O term in (!20|) which depends both on k and on r, and the fact that it becomes 
difficult for larger k, even with our cubic accelerant, to take sufficiently many primes 
p in the multivariate Euler product A^. The bound (128|) provides a rough estimate 
for the relative error introduced by truncating the Euler product, using our cubic 
accelerant at a given P, and very small e (strictly speaking, the remainder term 
in fl28l) needs to have the replaced by p-3+8fce take into account the choice of 
ej). Therefore, P 10"^'^^**^/^ is a crude estimate for how many primes are needed 
in the Euler product to achieve a relative error smaller than 10~^'S'*^. 

Rather than working with explicit constants in the truncation bounds, both for 
the O term in fl2Ul) and the Euler product, we experimented by taking different values 
of e = 10"°'^^*'^ and P, using our estimates as guides. We inspected the stability of 
our numerical values of Cr{k) by comparing those computed for a given P against 
those with P replaced by the first prime smaller than P/3, and only outputting the 
digits that agreed. It seems, from our tables, that the coefficients Cr{k) with mid- 
range values of r are more stable and converge faster with respect to P, especially 
for larger values of k. We did not explore the reason for this, but presumably the 
lower terms, beyond those resulting from our cubic accelerant, have comparatively 
smaller coefficients for those values of r for which Cr{k) converges faster. 

Numerical values of the coefficients for 4 < A; < 13 thus obtained are presented 
in the Tables [2]- El below in scientific 'e' notation, for example 1.2e — 3 = 1.2 x 10~^. 
High precision values of Cr{k) for k = 1,2, 3, can be obtained from [CFKRSj . 

It is also worth mentioning that all the digits of the coefficients computed in this 
manner agree, except in a few instances where the last decimal place differs slightly, 
with the results of the first method of [CFKRS2] (see Section [2]). That method has 
the advantage of producing high precision values of the coefficients, but is limited to 
relatively small values of r. We reran the program used in |CFKRS2] for k < 13 and 
r < 7 and display those values in Table [10] for comparison. 



13 



14 



T 


Cr(,4j 


Cr[b) 




U 




1 y( 1 i^nni non^^ooToi onnncCrt oa 
i.4iDUUiUzUDzz/oizUUyooe-z4 


o.izy4/o4uyi4yiyiize-4U 


i 


5.45(J14(J57oll71oo55yoooU5oe-ll 


'7 oon /( 1 o^c/^ /I n /I /I CI on cn/^o ^ oo 

7.ooU4iz75o4y445ioU5yooe-zz 


o.o(Jo7ozo(jyyzo444z4oe-o7 


o 
Z 


oo ^^on^^ o /I 1 ono 1 1 oo /I on nn 
o.zo/ /zyoo4/yizUoiioo4oy /e-uy 


1 ^^n^^n^^ oo 1 n^^ txonn onn /I ^ in 

i.i ^ yozooiyoozyuoouy4e-iy 


z.oU /yzU ( ( 114ooolzoyoe-o4 


o 
O 


o CiC^ 1 /101 'Tnnnon'yn /iixni^mo^ nv 

z.yo4ii4di ^yyyoy/y4oyDyioe-ur 
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1.603303769e+02 
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-3.78219665e+01 



Table 2: Coefficients Cr{k) for k — 4, 5, 6. For A; = 4, 5 we used e — 10 and 
P — 942939827, 180343651 respectively. For k — 6 we combined two data sets using 
e = 10-2^ P = 25501199 for < r < 29, and e = 10"^^ P = 608121859 for 
30 < r < 36. 
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Tabl6 3: Co6ffici6nts ior k — 7 truncating at P = 11015647, and using e = 10 
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Tabl6 4: Co6ffici6nts for A; = 8 truncating at P = 1212569, and using e = 10"^^ 



17 



T' 




/ 




T' 




n 

u 




1 
± 


fins74'^87'^p 1 1 n 


9 


S n'i19Qfi979p in? 


Q 
O 


1 T 79'^fi94nfip 1 (T? 




1 9^'^c\c\'\R7f,Qp inn 

l.ZOOUUOO ( Uc/c-lUU 


o 


1 n481497'?7fip Q7 


u 


7 1 A^fi'^l nn'^9o 
/ . i^ouoiuuozL-yo 


7 


4 n891 744^Qfip Q9 


Q 

O 


1 QQ41 Q9^9Qnp 


Q 


8 4f;Q49nf>n88p 87 


1 


1 fi'^SI 441 7'?P 84 

O.lOOOl^^l ( Oc-O^ 


1 1 
X X 


1 n4'^'^78QQ49p 81 


1 9 


O.UO / OlOUU^Oc/C / cJ 


1 


R 99414198Qfi1p 77 


1 4 
x^ 


1 QS'^l 487^71 fip 74 




A '^4Qnfi'iQ'^fi'^np-79 


1 R 


8 7n843439nQ7p-7n 


1 7 
1 1 


1 "10781 8788 'in9p-fi7 


lO 


9 fiQ'?9(SzLnnfi^1 1 p fi^ 

^ .Uc/O^UTiUUUtJ 1 1 1> UtJ 


1 Q 


4 1 898^^1 9Q479fifip fi'^ 


9n 

zu 


^ qq81 n841 '?^49fip fil 


91 




99 


Q 7S1 fi'^S'^fin7'^'^1 p ^7 




1 1 1 ^Qnfi'^fi^9894p ^4 


94 


1 1 S9SQ'^fiQQfi9Q9p ^9 


9^ 
zo 


1 1 fifi'^fi9n^77fi'^9np 


9fi 


1 n7n74Q4871 81 9p 48 


97 


Q 1 ^SSfi'^nQfil 7qQp 47 


98 
zo 


7 '?n41 n^Q447'^1 47p 4^ 


9Q 

Lo 


^ 4'^'^^9n^841 q7q4p A'\ 


sn 


771 8n'^n8Qfi9'^1 1 p-41 


tj± 




39 

OZ 


1 478984Q9081749P-37 




8 '^481 41 1 5859835p-36 


34 


4 4ni n701 35763Q61 p-34 


35 


9 1 657003966q394p-39 

Z. XUU 1 UUOZUUiyOZTiL. OZ 


ou 


Q Q44Q'^'^Q(S'^'^Q44(^p "^1 


^^7 


4 9(S008480Q'^89'^^p 9Q 


oo 


1 7ni ^8fi^q7q^984qnp 97 


oy 


U.OOOc/ZrOOc/c/ 1 Oc/OUc-ZrU 


4n 


9 1 Q^81 041 (SOO'^'^On 94 


41 

^ X 


7 n849fifi7^1 8qn7'^p 9'\ 


49 


9 19^9197809090740^ 91 


4'^ 


Q9941Q887098mp 90 


44 


1 ^n98^'^n(S4^i p i8 

L.OO lOUZOOOUU^O IC lO 




fi7nfi9n7fi49fifi784p-17 


4fi 




47 


1 fifiQ7'^1 ^1 44Q71 np-1 4 




3 1 5Q48978435nnQp-1 3 


4Q 


^n744Q8691 31 347p-1 9 


5n 


8 89746n86q45q55p-1 1 


51 


1.29834184462004e-09 


52 


1.74846308044295e-08 


53 


2.1508513613974C-07 


54 


2.4107241046116C-06 


55 


2.4552071751795C-05 


56 


2.265572158192C-04 


57 


1.888394621727e-03 


58 


1.4172609520338e-02 


59 


9.546112504952e-02 


60 


5.751640191248e-01 


61 


3.08994036050e+00 


62 


1.47569216973e+01 


63 


6.2480977623e+01 


64 


2.3393469844e+02 


65 


7.721276404e+02 


66 


2.23430645e+03 


67 


5.60365351e+03 


68 


1.1907448e+04 


69 


2.062257e+04 


70 


2.776900C+04 


71 


3.06185C+04 


72 


4.7187e+04 


73 


1.2202C+05 


74 


2.314e+05 


75 


1.255e+05 


76 


-4.658e+05 


77 


-1.07e+06 


78 


-5.794e+05 


79 


6.75e+05 


80 


8.27e+05 


81 


1.3e+05 











Table 5: Coefficients for A; = 9 truncating at P = 170741, and using e = 10 
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Table 6: Coefficients for A; = 10 truncating at P = 675929, and using e = 10 
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Tabl6 7: Co6ffici6nts for /c = 11 truncating at P = 85889, and using e = 10 
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Table 8: Coefficients for A; = 12 truncating at P = 12979, and using e = 10-^^ 
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Table 9: Coefficients for A; = 13 truncating at P = 1699, and using e = 10 
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Tabl6 10: For comparison, high pr6cision valu6S of Cr{k), A < k < 13, r < 7, 
comput6d using th6 program for m6thod 1 of [CFKRS] . Th6 valu6S h6r6 agr66 with 
thos6 in Tabl6s|2]-in]comput6d using our cubic acc6l6rant m6thod with 6xp6rim6ntally 
d6t6rmin6d r6maind6rs, 6xc6pt th6y ar6 occasionally slightly off in th6 last d6cimal 
plac6. 



5 Checking the Moment Polynomial Conjectures 



Let ^ 

Datafc(T)=/" |C(l/2 + 2t)|2^c^t, (32) 

with To = 14.134725 . . . being the imaginary part of the first non-trivial zero of zeta, 
and let 

Conjecturefe(T) = / Pk{\og—)dt. (33) 
Jto 27r 

We used our tables of Cr{k) and a simple integration by parts to compute the latter 
for given T. Note that because the leading coefficients of Pk{x) are very small, in the 
range of T considered, the Cr(fc)'s with midrange and higher values of r contribute the 
dominant amount to the integral. To accurately compute the prediction of |CFKRS] 
one does need the lower terms of the polynomial Pk{x). Table [11] displays, for 
1 < A; < 13 and T = 10*, the r for which the corresponding term in fl55]) contributes 
the dominant amount to integral. 



k 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


r 








2 


6 


11 


18 


26 


37 


49 


64 


80 


99 


119 



Table 11: The value of r for each k for which the corresponding term in ( l33l) con- 
tributes the dominant amount when T = 10*. 



To calculate we used the tanh-sinh quadrature scheme |Bal IBLJ] to accu- 
rately estimate each integral between consecutive non-trivial zeros of the zeta func- 
tion on the critical line. All our computations of f l32|) were carried out using 64 bit 
machine doubles. To tabulate all the zeros up to, and slightly beyond, T = 10*, we 
used Rubinstein's C++ L-function package Icalc [R]. It applies the Riemann Siegel 
formula to evaluate ({1/2 + it) and look for sign changes of the Hardy Z-function, 
Brent's method to compute the zeros of zeta |Brj . and a variant of Turing's test [E] 
to confirm that all zeros up to given height have been found. 

Figures [l]l2] depict the relative remainder term for 1000 values of T between and 
10*, roughly spaced apart every 10^. More specifically, we let Tj be the imaginary 
part of the first zero above lO^j, so that Tj ^ lO^j, and plot the values of 

Datafc(Tj) - Conjecturefc(Tj) 
Conjecture^ (Tj) ' 

for 1 < j < 1000, and 1 < A; < 12. 
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We also calculated the running average of the remainder term divided by the 
running average of the conjecture: 



Ei=i(Datafc(Tj) - Conjecturefc(T,)) 



■'J 

If we define 



Ej=i Conjecturefc(T,) 
SmoothDatafc(T) = /" \C{l/2 + it)\^''{l - t/T)dt, (36) 

JTn 



and 



t 

SmoothConjecturefc(T) = / Pk{\og—){l-t/T)dt, (37) 

Jto 27r 

then fl35|l gives a discrete approximation to the smoothed relative remainder: 

SmoothDatafc(Tj) — SmoothConjecture;j(Tj) 
SmoothConjecture;. (Tj) 

The reason for considering smoothed moments is that the noisy remainder terms of 
the sharply truncate moments, when averaged, tend to be smaller. 

Note that the vertical axes in these figures change from plot to plot to allow us 
to meaningfully display the relative remainder terms, which as a whole get larger, as 
k increases. We also set the zoom level to show the running averages in some detail. 
As a compromise, a few outliers are omitted from these plots for smaller T, roughly 
up to 10^, and k < 4. 

Table [12] lists the values of Datafc(T), and Conjecture;j(T) for = 1, . . . , 13 and 
T = 100000000.64, the first zero of zeta above 10^ We also list the values of the 
averages over all 1000 values of Tj 

1000 

^5^Data,(T,) (39) 
i=i 

and 

T 1000 

5^Conjecturefc(T,). (40) 



1000 

j=i 

Our data supports the CFKRS conjecture for the full asymptotics of the moments 
of zeta as described in equation ([T]), though for larger k, it is difficult to gauge the 
size of the remainder term. 
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For k = 1 the data suggests an even stronger remainder term of O (T^^'^'^'^) , 
supported by the agreement of Datai(T) with Conjecturei(T) to about 3/4 of the 
decimal places left of the decimal point for the values of T examined. The relative 
remainder term is of size around 10~^ when T ^ 10^, and, typically, a couple of 
orders of magnitude smaller when averaged. 

For k = 2 the data agrees with the conjecture to about half the decimal places, 
with a relative remainder term of size around 10^^. For k = 3, the agreement is to 
slightly less than half the decimal places. 

For fixed T, as k increases, the moments have the effect of amplifying the largest 
values of |C(l/2 + 'it)\. This can be seen in our plots, for larger k, where the re- 
mainder terms are qualitatively the same, with large jumps at the same values of T 
corresponding to relatively large values of zeta. It therefore becomes more difficult 
to ascertain, as k grows, whether an upper bound for the remainder term of the 
form Ofe(T^/^+'^) holds. Nonetheless, the table and figures reveal an excellent fit for 
the CFKRS prediction with the moments which persists through to the 24th and 
26th moments, where the relative agreement is to within around one to two decimal 
places. 

In some sense, the fit between columns 2 and 3 of Table [12] is better than it 
ought to be for larger k, for example more than three decimal places for k = 13, but 
only agreeing to around 90% for the average remainder. A quick inspection of the 
figures reveals that the relative remainder, at T = 10®, happens to, fortuitously, best 
its neighbours, especially for larger k. Nonetheless, the overall agreement between 
the CFKRS prediction and our data across all values of k and T, as depicted in 
the figures, lends strong support to their conjecture for the full asymptotics of the 
moments of the zeta function. 
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3.43851173e+32 


3.43090327e+32 


7.0665341e+30 


7.13020297e+30 
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1.72388578e+36 


1.71918466e+36 


3.26838252e+34 


3.32617445e+34 


10 


9.2785049e+39 


9.251733e+39 


1.6228859e+38 


1.6729912e+38 


11 


5.2991086e+43 


5.2863072e+43 


8.5447722e+41 


8.967058e+41 
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3.182548e+47 


3.179455e+47 


4.726347e+45 


5.076008e+45 


13 


1.995625e+51 


1.999377e+51 


2.726982e+49 


3.013406e+49 



Table 12: The values of Datafc(T) and Conjecturefc(T) at T = 100000000.64, the first 
zero after 10®. We also display, in the last two columns, the values of the averages 
over all 1000 values of Tf Yl!j^=i Datajt(Tj), and Ylf=i Conjecture;j(Tj). 
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Figure 1: Plot of the relative remainder f l34p . and running average of the remainder 
(1551) . for 1000 values of T, and 1 < A; < 6. Horizontal axis is T. 
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Figure 2: Plot of the relative remainder f l34p . and running average of the remainder 
for 1000 values of T, and 7 < k < 12. Horizontal axis is T. 



